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TECHNICAL NOTE 5^25 

SOME CALCULATIONS OF THE LATERAL RESPONSE OP TWO 
AIRPLANES TO ATMOSPHERIC TURBULENCE WITH RELATION TO THE 
LATERAL SNAKING IROBLEM^ 

By John D. Bird 

SUMMARY 


Calculations are made of the lateral response to representative 
time histories of atmospheric tiorhulence for two airplanes having 
widely different dynamic properties, and explanations for their differ- 
ences in behavior are given. 

The resiilts of the calculations indicate that, under the proper con- 
ditions, atmospheric turbulence can initiate and maintain a lateral 
hunting oscillation of an airplane, and that this oscillation can be 
fairly regular in both amplitude and frequency. This effect is more 
pronounced for lightly damped airplanes. It is felt that this phenome- 
non may be the cause for some of the cases of airplane snaking that 
have not been explained by other considerations. 


INTRODUCTION 


Several recent hi^-speed airplanes exhibit a tendency to develop 
and maintain lateral hunting oscillations of roughly constant amplitude, 
which are generally referred to as snaking oscillations . Tbese oscil - 
lations are of essanti aily t he same neriod as the normal Dutch roll 
oscillat ion nf ■ the . ai mlan a,^d .usually have ni npHtud es of the order 
of '^1'^ yaw . They are considered to be a nuisance to the pilot in that 
they cause loss of confidence in the airplane's response to control and 
make the airplane less satisfactory as a gun platform, 

Adeqviate explanations have been offered for this behavior in 
specific cases; however, there are still numeroTis occvirrences for which 
no satisfactory explanation has been advanced. Some of the explanations 
for this motion are associated with nonlinear aerodynamic character- 
istics which result in different rates of damping for the large and 

T — — - - - 

Supersedes the recently declassified NACA RM L50F26a, 1950. 
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small amplitude-ranges of motion. An example of this nature could 
arise frcm a poor fairing at the Juncture of the tail surfaces. Other 
causes of snaking can be associated with small amovints of slack in the 
rudder control system, or the effects of fuel sloshing. 

Some impublished experiments conducted In the Langley stability 
tunnel have suggested the possibility that motionSj^ similar to those 
described as snaking, may result from the turbulence which exists in 
the atmosphere. The present paper, therefore, constitutes a preliminary 
investigation of this possibility. Calculations are made of the lateral 
response to representative time histories of atmospheric turbulence for 
two airplanes, having widely different dynamic prd^rties, and exjilana- 
tlons for their differences in behavior are given. " Calculations suad 
experimental records of the response of a model, which has freedom only 
in yaw, to the roughness in the stability- tunnel air stream are given 
for illustrative purposes. 


SYMBOLS 


The stability system of axes is used ln_ the" present analysis. This 
axis system has its origin at the center of gravity. The Z-axis is in 
the plane of“ symmetry and perpendicular to the relative wlnd^ the— X-axis 
is in the plane of symmetry and perpendicular to the Z-axis, and the 
Y-axis is perpendicular to the plane of symmetry. The positive direc- 
tions of the stability axes and of angular displac^ents of the airplane 
are shown in figure 1. The coefficients of forces and moments employed 
in this paper are in standard NACA form. 

The coefficients and symbols are defined"as“follows: 

Cp lift coefficient, Lift/qS 

Cy lateral -force coefficient, Y/qS 

C^ rolling -moment coefficienty L/qSb 

Cn yawing -moment coefficient, N/qSb 
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^XZ 

^XZ 

^z 
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p 




radius of gyration about X-axis, 
radius of gyration about Z-axis, 
product-of -inertia factor, ’]Sy^/\P- 
radius of gyration about X-axis 
radius of gyration about Z-axis 
product-of -inertia factor 
relative density factor, 
moment of inertia about Z-axis 

angle of bank of airplane, radians unless otherwise noted 
azimuth an^e of airplane, radians unless otherwise noted 
angle of sideslip of airplane, radians unless otherwise 
noted, tan”^ ~ 


°2V 

°2V 

^rb 


dC 


n 


N rb 
°2V 


dimensionless 


dlmens ionles s 


dimens ionles s 



h 
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CT 

o’o 

lOCX) 

f 

t 

Db 



Ti inclination of principal longitudinal axiS“Of inertia with 

respect to flight path} positive when the principal axis 
is above flight path at nose 

W weight 

S wing area 

h wing span 

A aspeci~ratio 

g acceleration due "to gravity 

X longitudinal force along X-axis 

Y lateral force along Y-axis 

Z normal force along Z-axis, Lift = - Z 

L rolling moment about X-axis 


angle of air stream with- respect to initial fll^t-path 
direction of airplane, radians unless otherwise noted 

amplitude of oscillation in air-stream direction 
amplitude of oscillation in model heading 

relative an^jlification, V in feet per second 

frequency, cycles per unit time 
time ^ 

differential operator, d/ds^ 
dimensionless time, tv/b 
a partlcTjlar time 


N 


yawing moment about Z-axls 
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pb/ 2 V wing-tip helix angle, radians 

rb/2V yawing -velocity parameter, radian measure 

p rolling angular velocity about X-axis, radian measure 

r yawing angular velocity about Z-axis, radian measure 

V linear velocity of airplane along Y-axis 

Vg^g^ lateral gust velocity with respect to undisturbed position 
of airplane 

V free -stream velocity 

q. dynamic pressxire, ^ 

p mass density of air 

P period of free lateral oscillation 

T 1/2 time for free lateral oscillation to damp to one-half 

amplitude 


CALCULATION METHODS 


General methods are given in references 1 and 2 for calculating 
the lateral response of airplanes to gusts. Reference 1 indicates that 
the response of an airplane to an arbitrary gust structure may be 
obtained by superposition of solutions for unit gusts, which, in the 
limiting case of a continuous disturbance function. Involves the evalua- 
tion of Duhamel’s integral. It is also pointed out in reference 1 that 
a rigorous analysis of gust effects requires consideration of penetration 
time and of the aerodynamic lag in building up the lift on the surfaces. 
Examination of some of the penetration effects indicated that their 
magnitudes, were Fnnn .1 1 compared with the effects of sideslip. Eie rolling 
component of .the measured turbiilence considered for this paper was, of 
course, unlmown; and, thus, all gust disturbances were necessarily con- 
sidered to be in a single plane. For the purpose of the calculations 
of this paper, the turbulence was assumed to contribute nothing more than 
an effective change in sideslip of the airplane. Thus, in a side gust 

dC^^ dC 7 

of velocity Vcr, the angular and linear disturbances are a — a — t, 

do dcr 
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dC, 


and 0 - — -wixere the derivatives are numerically equal to CL , > 

do P P- 

dCy 

and Cy , respectively. In this analysis a — ^ was assumed to he zero, 
■^p d0 


The lateral response of two airplanes to representative time his- 
tories of atmospheric turbulence was calculated by obtaining the motion 
of the airplanes following the application of unil;’ yawing- and rolling- 
moment coefficients and then by evaluating Duhamel‘s integral with this 
unit solution as the response variable and the record of the lateral 
fluctuation in air-stream direction as the forcing function. 


Now Duhamel * s integral may be written 







where - s-j^^ is the lateral response of the airplane to a 

unit a. This equation may be broken into two Integrals: 


where >lfjj and - s-j^^ are the lateral responses ofi 

the airplane to unit yawing- and rolling -moment coefficients, respec- 
tively, and dCjQ^/d 0 and dC^/d0 are numerically equal to Cj^ 


and C^p, respectively, 
as a single integral: 


The last equation, of course, may be expressed 


^s-b3_ 






^(sbi - St) 



The solutions to the lateral equations of motion following the applica- 
tion of unit yawing- and rolling -moment coefficients were obtained by 
use of an automatic digital computing machine and the procedures of- 
reference 2. Duhamel *s Integral was evaluated by a numerical integra- 
tion in a manner similar to that of ref^ence 5 to obtain the motion 
.of the airplane in. response to the turbiilence. This calculation was 
also carried out on an automatic computing machine. 

A brief resume of the methods used for calculating the lateral 
frequency -re spouse characteristics of the airplanes considered is given 
In the appendix along with the equations of motldh from which the 
response of the airplanes to unit disturbances was calculated. 
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RESULTS ARD DISCUSSION 
CalcTjlated Airplane Motions in Turbiilent Air 


The mass and aerodynamic characteristics of two airplanes having 
widely different operating conditions and dynamic characteristics are 
given in table I. Airplane A is a low-speed, low-altitude airplane j and, 
airplane B is a high-speed, high-altitude research airplane, which is 
known to exhibit small continuous lateral oscillations imder certain 
flight conditions . 

Figures 2 and 5 show the results of the calculation of the lateral 
response of airplanes A and B to two known distributions of atmospheric 
turbulence. Ihe turbulence shown in figure 3 "was measured by means of 
a sensitive recording, pitot-static tube mounted on an airplane which 
traversed a region of turbxiLent air and so recorded the fluctuations in 
forward speed through the region. The txarbulence shown in figure 2 was 
measured by means of an accelerometer mounted at the center of gravity 
of an airplane which traversed a region of turbulent air. Both of these 
records of turbulence were considered to be fluctuations in sidewise 
velocity for these calculations. Reference 4 gives information on the 
measurement of almiospheric turbulence and justification for the assump- 
tion that the turbulence is isotropic. 

The gust distributions shown In the figures were assumed to exist 
in like fashion along the flight paths of both airplanes. The high- 
speed airplane, of course, encounters gusts with a greater frequency 
than the low-speed airplane. The gust velocities are assumed for these 
calculations to be the same for all altitudes. These results are thus 
of a qualitative nature. 

The distance traversed by the airplanes is used as an abscissa in 
the plots given, and the azimuth angles of the airplanes are chosen to 
indicate the oscillation perfoimed. The azimuth angle should be roughly 
proportional to the apparent lateral movement of the horizon. 

The motions of the two airplanes in response to the two regions of 
atmospheric turbulence are markedly different (figs. 2 and 3)- Air- 
plane A shows a response which mi^t logically be termed by the pilot 
as rough airj that is, the airplane responds in almost direct proportion 
to the local gustiness and subsides to little or no motion as the gusti- 
ness subsides. The response of most airplanes in the past seems to have 
been of this nature. For example, see the small amplitude motions meas- 
ured in flight during the investigation reported in reference 5 . Air- 
plane B shows a response which bxiilds up to a fairly steady lateral oscil- 
lation, which is almost Independent of the local txirbulence. This motion 
is very similar in character to motions which have been termed snaking. 
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The-periods of the lateral oscillations, shown in figures 2 and 3, 
are close to the period of the classical free-lateral oscillations of 
the airplanes (table I) . This is more nearly the case for airplane B_ 
than airplane- A. 

Figure 4 shows the motions of both airplanes compared with the 
angular motion of the air with respect to the undisturbed attitude of- 
the airplane for the case of the turbulence determined from airspeed 
fluctuations. It should be noted that the turbulence record corresponds 
to different amplitudes of a for the two airplanes because of differ- 
ences in forward speed. It is easy to note a closer approximation in 
the case of airplane A to a one-to-one correspondence of air motion to 
airplane motion than for the case of airplane B, where there is little 
apparent relation between the air and airplane motions. It should be 
mentioned here that the long period change in heading shown for air- 
plane A in figure 3 was modified to some degree in the preparation of 
figure 4^ in order to have the air and airplane motions oscillate about 
the same mean and so make for an easier comparison of- the two motions. 

A comparison of a part- of the lateral motion of airplane B, as 
calculated from the atmospheric turbulence (fig. with the snaking 
of this airplane recorded during a flight test is shown in figure 5* 

This comparison merely confirms the statement made previously that 
lateral motions arising from turbulence in the air can be similar in 
nature to flight measurements of snaking motions. The flight conditions 
for the two motions given are— not-dLdentical, and the atmospheric turbu- 
lence existing during the fli^t test is unknown. The indications are, 
however, that turbulence having about one -third the magnitude of that~ 
employed for figure 3 W 9 uld be required to maintain a hunting motion of- 
alrplane B comparable with the snaking motion that it exhibited in fU^t. 
It should be pointed out tliat the disturbances do not have to be of the 
type generally referred to as sheirp-edged gusts as shown in figure 2, 
but may be of a more gentle nature, as shown in figure 3. 


Free Motions of Ifodel in Wind Tunnel 

As an example of the response of a free body to turbiilence in the 
air stream, figure 6 gives the experimental and calculated hunting 
motion of a model mounted with freedom only in yaw in the air stream of 
the Langley stability-tunnel test section. The measured time history 
of the air -stream azimuth angle from which the model motion was calcu- 
lated is also shown. The calculation procedure was much the same as has 
been given previously for airplanes A and B, The time history of-alr- 
stream direction was obtained by use of-a recording electronic pitot emd, 
although the percentage error in the magnitude of the air-stream angles 
may be fairly large, the nature of the fluctuations should he accurate. 
The e3q>erimental model motion was not recorded at the same time as the 
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air-stream azimuth angle. This fact does not change the nature of the 
result, however, because the air-stream turbulence was of a similar 
nature for a number of different recordings. 

The tunnel model was free only in yaw, was mounted on flexure plate 
in order to minimize friction, and consisted of oaxLy a fuselage and 
vertical tail. The mass and dimensional characteristics of the model 
oscillating system are given in table II. 

It can be seen that the experimental and calculated motions of 
figure 6 are of a similar nature, indicating that the air-stream turbu- 
lence is a significant factor in the hunting motion experienced by the 
model. 


Frequency -Response Characteristics 

The response of airplanes A and B, and of the tunnel model mounted 
with freedom in yaw, to sinusoidal forcing functions of various fre- 
quencies is given in figure 7 * Hie forcing functions are in the form 
of changes in heading of the approaching air stream. The results are 
given in terms of the amplitude of motion in radians induced by a sinus- 
oidal lateral gust distribution having an amplitude of 1000 feet per 
second and in terms of the phase lag of the motion behind the forcing 
function. 

An examination of these c^^rveB indicates the source of the differ- 
ence in response of airplanes A and B to atmospheric tiorbulence. Air- 
plane B has a frequency response not lonlike the characteristics of an 
electronic band-pass filter which excludes those harmonics of the 
applied frequency which are very much different from the resonant fre- 
quency. A good measure of the selectivity of response of these airplane 
is the ratio of the amplification at the natxaral frequency to the ampli- 
fication at very low frequencies. Airplane A responds to a greater 
degree than airplane B to those frequencies that are different from the 
resonant frequencyj thus, the tendency for the one-to-one correspondence 
of air direction to airplane azimuth angle shown for airplane A in fig- 
ure 4 . In general, the sharper the frequency-response curve the more 
nearly the response to atmospheric turbulence approaches a sinusoidal 
motion. 

Any effect which reduces the rate of free damping of an airplane 
should tend to increase the peak of the frequency-response curve and 
make the phase-angle shift at the natural frequency more abrupt. Air- 
plane B, of coiirse, has a low rate of dancing for the condition inves- 
tigated herein (table I). 

The rate of free damping can be affected to a marked degree by a 
change in one or more of the aerodynamic stability derivatives. This 
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fact voxild make calcxolations, ‘based on estimated derivatives, of a 
questionable nature unless some form of check is available. The sta- 
bility derivatives used for the calculations of this paper are believed 
to be reasonably accurate because the experimental and calculated rates 
of ^ lateral oscillation compare veil (table I). 

The frequency -response curve of the oscillating model mounted in 
the air stream of the Langley stability tunnel is between the curves of 
airplanes A and B with regard to select ivi-ty but is large compared with 
both as regards over -all response (fig. 7) . 


CONCLUDING RIMAEKS 


The results of the calculations of the lateral response of two air- 
planes to atmospheric turbu3.ence indicated that, under the proper con- 
ditions, atmospheric- turbulence can initiate and maintain a lateral 
h-unting oscillation of an airplane, and that this oscillation can be 
fairly regular in both amplltiade and frequency. This effect- is more 
pronounced, for lightly damped airplanes. It is suggested that this 
phenomenon may be the catzse for same of the cases of airplane snaking 
that have not been es^lained by other considerations. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Yo.,, June 2,'J, 1950* 
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APPENDIX 

CALCULATION OF LATERAL FREQUENCY RESPONSE 


The lateral frequency response of an airplane to an imposed sinus- 
oidal veirlation of wind direction may be calculated by solving the 
standard lateral equations of motion of an airplane (with the proper 
forcing terms added) for the steady-state motion. For a unit sinusoidal 
variation in air-stream direction the amplitude of this motion becomes 
the amplification factor or the amount that the air-stream fluctuation 
is magnified. 

Within the limits of the approximations discussed^ the lateral 
equations of motion of an airplane experiencing a sinusoidal variation 
in air -stream direction are 

/ 221 \ I 21 ^ fb 

^ “2 l^p^^xzDb - 2 == T 

■ I ^np^)^ + ^ 2jt ^ s-jj 

- C^cp + 2pbl>b^ + ^2|u^D^3 - Cyp) 3 = 0 

where the terms Cxp and Cnp on the right side of the first two equa- 
tions are considered to be the equivalents of and 

dq dq 

The variation of side force with rolling and yawing velocity, a 
term associated with the glide-path angle, and the side-force forcing 
function are omitted in these equations. Calcxilations showed these 
factors to be of little importance. Replacing the two right-hand terms 
of these equations by unity gives the equations from which the unit 
solutions were obtained for use in the calculations of response to 
arbitrary turbulence. 

Solving the equations for the steady-state motion gives for the 
az imuth angle i|/ 


= \|a^ + b^ sin 
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where € is the angle of lag of t he motion behind the disturbance , and 

for cTq = 1.0, the term \/a^ + b^ is the amplification factor for the 
given imposed frequency or ^ . 

= \/a2 + t.2 

Now the amplitude of the yawing motion for a given amplitude of lateral 
gust velocity is 


'I'o = 



This term is called the relative amplification when 1000 is substituted 
for and V Is given in feet per second. Now 


- ^u^u 



■ and 


where 


V _ t k^Ce 

Cu^ Ce^ 

^e = 

= mgfgB + mifg 

Ce - - Wb® 

Cu - - 


and 




2jtfb 

V 


mo = 0 
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”1 = (I °Zp%°np - I \%%) 

“2 = ■*■ - (l^fe'^n.p + ^HfcKxzPYgjoj^Oo 

“3 “(■ ■*■ 

Th.e expressions for the coefficients of the lateral-stability equa- 
tion A, B, C, D, and E are given in reference 6. 


The lag angle is given by 

-1 a 


e = cos 


\/a2 + b^ 


= sin' 


-1 


Va^+ b^ 


Substitution of various values of the imposed frequency in the previ- 
oiisly given expressions gives the frequency-response ctirve. 
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TABLE I 

MASS AKD AERODYNAMIC CHARACTERISTICS OF AIRPLANES 








^XZ 

"b, feet . . . 
Weight, pounds 
S, square feet 

• • • V ' . • 


lip 





* . . i . 



cyp 

V, feet per second . . . . 

Altitude, feet 

P (calculated) , seconds 
P (fligfht), seconds . . . 
T1/2 (calculated), seconds 

T1/2 (flight), seconds . . 


Airplane A 
16.8 
0.0061 

0.0264 

0.0058 

33.6 

8,700 

250 

0.551 

-0.280 

-0.085 

0.090 

-0.270 

- 0.049 

0.097 

-0.665 

257 

7,500 

3.7 

3.5 

2.5 

3.0 


Airplane B 

106.3 

0.0051 

0.0409 

-0.0006 

28 

11,050 

130 

0.343 

-0.474 

0 

0.224 

-0.170 

- 0.101 

0.217 

-0.878 

746 

30,000 

1.48 

1.80 

6.56 

6.50 
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TABLE II 

MASS AND AERODYNAMIC CHARACTERISTICS OE MODEL MOUNTED 
WITH FREEDOM IN YAW IN AIR STREAM OF THE 
LANGLEY STABILITY TUNNEL 


I 2 , foot-pouna-seoond2 

poTmd per square foot 

f feet 

V, feet per second . . 

*S^ square feet .... 
Cj^p per radiaji .... 


P (calculated), seconds . 
Ti/o (calcvilated) , seconds 


0.3 

itO 

3 

186 

2.25 

0.063 


-0.13 

0.8 

k.k 


*The symbols b and S are given as dimensions upon idaich 
aerodynamic coefficients are based for this model. 



Relative ulnd, V 


T 

z 

Section A-A 



Figure 1.- Stalsillty system of axes. Positive values of forces, moments 
velocities, and angles are indicated ty arrows. 
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Figure 3 *- Calculated lateral response of airplanes A and B to atmospheric 
turbulence . Turbulence obtained from eiirspeed fluctuation records. 
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Figure 4.- Ccmrparlson of lateral motions of airplanes A end B vith 
fluctuation in air-stream direction. Turljulence dietribitlon of 
figure 3. 
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(b) Calculated response of airplane to atmospheric 
turbulence (from figure 3) , 

Mach no.^0.74; altitude ,30^000 feet ^ weighty 1 1^000 pounds , 

Figure 5.- Comparison of calculated hunting motion of airplane B In 
response to atmospheric turbulence vith measured hunting motion of 
airplane B obtained from a flight test. 
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Figure 6.- Ccmparison of huntiiig motion of model mounted with freedom 
in yaw in at ability- tunnel air stream with calculated hunting motion 
determined from a representative ^ record of air-stream turbulence, 
q s: liO pounds per square foot. 
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(a) Relative amplification. 


of free oscillation 



Figure 7.- Lateral- frequency-response characteristics of airplanes A 
and B and model mounted with one degree of freedom in stability- 
tunnel air stream. 


NACA TN 



XACA-Lugley - B-IT-H - 1000 


Airplane B 
Airplane A 

Model in funnel^ freedom m yaw 


ro 

-p- 



Ratio of impressed frequency to natural frequency of lateral oscillation 
(b) Phase lag of motion behind impressed moment 

Figure 7-- Concluded, 


BtL VDVH 



